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Abstract 

In view of the recent start of the NA62 experiment at CERN that is expected to 
measure the —>■ tt+i/P branching ratio with a precision of 10 %, we summarise the 

present status of this promising decay within the Standard Model (SM). We do like¬ 
wise for the closely related Kl —t which will be measured by the KOTO exper¬ 

iment around 2020. As the perturbative QCD and electroweak corrections in both 
decays are under full control, the dominant uncertainties within the SM presently 
originate from the CKM parameters \Vcb\, |I4fe| and 7 . We show this dependence 
with the help of analytic expressions as well as accurate interpolating formulae. Un¬ 
fortunately a clarification of the discrepancies between inclusive and exclusive deter¬ 
minations of 1 14 ;, I and 1 14 b I from tree-level decays will likely require results from the 
Belle II experiment available at the end of this decade. Thus we investigate whether 
higher precision on both branching ratios is achievable by determining | 14 b|, | 14 b| 
and 7 by means of other observables that are already precisely measured. In this 
context ek and AMg^d, together with the expected progress in QCD lattice calcu¬ 
lations will play a prominent role. We find B{K~^ —)• vr+z/i/) = (9.11 ± 0.72) x 10“^^ 
and B{Kl = (3.00 ± 0.30) x 10“^^, which is more precise than using aver¬ 

ages of the present tree-level values of lUcbl, |Uub| and 7 . Furthermore, we point out 
the correlation between B{K^ —?■ Tr^nn), B{Bs —)• pAgi~) and 7 within the SM, that 
is only very weakly dependent on other CKM parameters. Finally, we update the 
correlation of —)• Anv with the ratio s'/e in the SM taking the recent progress 

on e'/e from lattice QCD and the large N approach into account. 
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1 Introduction 


The measurements of the branching ratios of the two golden modes — )• ■k'^uv 

and Kl —)■ will be among the top highlights of flavour physics in the rest of 
this decade. K~^ —?■ is CP conserving while Kl —)• is governed by CP 

violation. Both decays are dominated in the SM and in many of its extensions by 
Z penguin diagrams. These decays are theoretically very clean, and the calculation 
of their branching ratios within the SM includes next-to-leading order (NLO) QCD 
corrections to the top quark contributions NNLO QCD corrections to the 

charm contribution [4]-[^ and NLO electroweak corrections [7[|^ to both top and 
charm contributions. Moreover, extensive calculations of isospin breaking effects 
and non-perturbative effects have been performed [^, 11 . Reviews of these two 
decays can be found in 12 -18 and their power in probing energy scales as high as 
several hundreds of TeV has been demonstrated in |19|. 


In view of the recent start of the NA62 experiment at CERN that is expected 
to measure the —)• vr+nn branching ratio with a precision of 10 % compared to 

and the expected measurement of Kl —)• by KOTO 

it is the right time to summarise the present 


20 


21 


the SM prediction 
around 2020 at J-PARC 1151 22 


status of these decays within the SM. This is motivated in particular by the fact 
that different estimates appear in the literature due to different inputs used for the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, which presently constitute 
the main uncertainty in the SM predictions for these two branching ratios. This has 
been stressed in |23| , where the dependence of both branching ratios on the chosen 
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values of |I46| and \ Vub\ extracted from tree-level decays has been studied (see Table 
3 of that paper). 

At this point two strategies for the determination of the contribution of the SM 
dynamics to these decays are envisaged; 


Strategy A: The CKM matrix is determined using tree-level measurements of 


IK. 


iKfel, iKbl, 7, 


( 1 . 1 ) 


where 7 is an angle of the unitarity triangle (UT). As New Physics (NP) seems to 
now be well separated from the electroweak scale, this determination of the CKM 
matrix is not expected to be polluted by NP contributions^ Inserting these inputs 
into the known expressions for the relevant branching ratios (see Section then 
allows a determination of the SM values for the K —)• ttuv branching ratios inde¬ 
pendently of whether NP is present at short distance scales or not. The departure 
of these predictions from future data would therefore allow us to discover whether 
NP contributes to these decays independently of whether it contributes to other 
decays or not. This information is clearly important for the selection of successful 
extensions of the SM through flavour-violating processes. 

Unfortunately, this strategy cannot be executed in a satisfactory manner at 
present due to the discrepancies between inclusive and exclusive determinations of 
iKbl and iKfel from tree-level decays. Moreover, the precision on 7 from tree-level 
decays is still unsatisfactory for this purpose. While the measurement of 7 should be 
significantly improved by LHCb in the coming years, discrepancies between inclusive 
and exclusive determinations of |Kb| and |Kfe| from tree-level decays are likely to 
be resolved only by the time of the Belle II experiment at SuperKEKB at the end 
of this decade. 

The clarification of the discrepancies between inclusive and exclusive determina¬ 
tions of iKfel and I Kb I from tree-level decays is important, but there are reasons to 
expect that the exclusive determinations will eventually be the ones to be favoured. 
First of all, exclusive measurements are easier to perform than the inclusive ones. 
Equally important, due to the significant improvement in the calculations of the 
relevant form factors by lattice QCD, exclusive determinations are more straight¬ 
forward than the inclusive ones. This is opposite to the philosophy of ten years ago, 
where QCD lattice calculations were still at the early stage and inclusive determi¬ 
nations were favoured. 

Yet, from the present perspective it is useful to study the SM predictions for 
K+ —)■ -K^vv and —)• in the full range of |Kb|) |Kb| and 7 known from 

tree-level decays, as this will clearly demonstrate the need for the reduction of 
parametric uncertainties. This will also allow the SM predictions for these decays 
to be monitored as the determination of 7 will improve in the coming years at the 
LHC. This should be of interest in view of the first results on K~^ —?• from 

NA62, which are expected already in 2016. Moreover, it will also be of interest to see 
how other observables, like ek, AM*, AM^, and rare Bg d decays are modified when 


TT ' vv 


the parameters in ( 1 . 1 ) are varied, and what their correlations with K~^ 
and Ki ^ are within the SM. As we will see, some of these correlations are 
practically independent of |Kb| and |Kb| and as such are particularly suited for a 
precise tests of the SM. 


^Recent analyses of the room left for NP in tree-level decays can be found in 24 -26 . 
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Strategy B: Here the assumption is made that the SM is the whole story and 
all available information from flavour-changing neutral current (FCNC) processes is 
used to determine the CKM matrix. Our strategy here will be to ignore tree-level 
determinations of \Vub\ and \Vcb\, as the discrepancies mentioned above could also 
result from experimental data, which will improve only at the end of this decade. 
Similarly, the tree-level determination of 7 will be left out. Then the observables to 
be used for the determination of the CKM parameters will be[^ 

ek, AMs, AMrf, S^Ks^ (1-2) 


accompanied by lattice QCD calculations of the relevant non-perturbative param¬ 
eters. In this manner also \Vcb\, Ifufel and 7 can be determined. This is basically 
what the UTfit 


27 and CKMfitter |28| collaborations do, except that we ignore 


the tree-level determinations of \Vcb\, \Vub\ and 7 for the reasons stated above. As 
the dominant top quark contribution to ek is proportional to \Vcb\^ and AMg^d 
are proportional to |I 4 bp, a useful determination of |I 4 fe| can be obtained from 
these quantities]^ The full UT is then constructed by using the ratio AMd/AMg 
and S^Ks- We find that with the most recent lattice QCD input on the parameter 
^ 1 ^, the determination of 7 in this manner is impressive, and also the value of \ Vcb\ 
is significantly more accurate than from tree-level decays. In the case of \Vub\ the 
accuracy is found to be comparable to the most recent exclusive determination [32] . 

It should be emphasised that while strategy A is ultimately the one to use 
to study extensions of the SM, the virtue of strategy B at present is the greater 
accuracy of the SM predictions for the observables that we consider. By simply 
imposing constraints from several measurements we arrive at narrow ranges for the 
parameters in (1.1) - given that the SM is the whole story. 


In the present paper we will follow these two strategies using the most recent 
inputs relevant for both of them, in particular the ones from lattice QCD. In Sec- 
tion[^we summarise the present status of the K~^ —>■ tt^vv and Kl —)> tt^vD decays 
in the SM and discuss the main uncertainties with the help of analytic expressions. 
In Sections and we follow strategies A and B, respectively, and present in some 
detail our numerical results. In Section we present an updated analysis of the 
correlation of —)• and the ratio e'/e in the SM. We conclude in Section]^ 

In the appendices we collect a number of additional expressions that we used in our 
analysis. 


2 Basic formulae 


We present here the basic formulae for the branching ratios for the K~^ —?■ 
and Kl —)• decays in the SM. This section can be considered as an update to 
the analogous section (Section 2) of 12 , a review of these decays from 2007. The 


main advances in the last eight years are: 

• computation of complete NLO electroweak corrections to the charm quark 
contribution to —>■ in |^; 


^Note that the present determination of has no impact on the CKM parameters in the SM. 

^The strategy for the determination of \Vcb\ from ek is not new 29 and has been considered recently 


in 30 . See in particular formula (29) in that paper. 
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• computation of complete NLO electroweak corrections to the top quark con¬ 
tribution to K~^ — 7 - and Ki —)• in [^; 

• reduction of uncertainties due to rndmc) and as{Mz), with the last 

two relevant in particular for the charm contribution to —)• 


While incorporating these advances in our presentation we will also include 

• NLO QCD corrections to the top quark contributions (H and NNLO QCD 
corrections to the charm contribution 


isospin breaking effects and non-perturbative effects 10,11 


2.1 K+ 




The branching ratio for K~^ —)■ in the SM is dominated by penguin di¬ 

agrams, with a significant contribution from box diagrams. Summing over three 
neutrino flavours, it can be written as follows [3,11 


B{K^ —^ ■K^vu) = K+(l -I- Aem)- 


ImA< 

A5 


X{xt) 


+ 


RsAc -p (v\ ! RsAt , , 

-^Pc{X) + -^X[xt) 


( 2 . 1 ) 


with 


K+ = (5.173 ±0.025) • 10 


-11 


A 


0.225 


Aem = -0.003. (2 .2) 


Here xt = A = \Vus\, K = V*,Vid are the CKM factors discussed below, 

and K+ summarises the remaining factors, in particular the relevant hadronic matrix 
elements that can be extracted from leading semi-leptonic decays of K'^, Ki and Kg 
mesons [^. Aem describes the electromagnetic radiative correction from photon 
exchanges. X{mt) and Pc{X) are the loop functions for the top and charm quark 
contributions, which are discussed below. An explicit derivation of (2.1) can be 
found in |^. The apparent large sensitivity of B{K^ —)• to A is spurious as 

Pc{X) A ^ (see (|2.6[)) and the dependence on A in \2.2\ cancels the one in (| 2 .l[) 


to a large extent. Therefore when changing A it is essential to keep track of all the 
A dependence. _ 

11 , the MS scheme with 


In obtaining the numerical values in (2.2) 


Wey,{Mz) = 0.23116, 


a{Mz) = 


1 


127.925’ 


(2.3) 


has been used. As their errors are below 0 . 1 % these can currently be neglected. 
Note, however, that although the prefactor of the effective Hamiltonian, a/ sin^ 9^, 
is precisely known in a particular renormalisation scheme (MS in this case) it re¬ 
mains a scheme dependent quantity, with the scheme dependence only removed by 
considering higher order electroweak effects in A —)• ttz/P. An analysis of such ef¬ 
fects in the large ruf limit [^ demonstrated that in principle this scheme dependence 
could introduce a ±5% correction in the K —)• vrz/P branching ratios, and that with 
the MS definition of sin^ 9w these higher order electroweak corrections are found 
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below 2%. However, only the complete analysis of two-loop electroweak contribu¬ 
tions to if —7- TTvv in for the top contribution could put such expectations on firm 
footing. The same applies to the NLO electroweak effects in the charm contribution 
to —>■ evaluated in [^. 

The short distance function X[xt) relevant for the top quark contribution, in¬ 
cluding NLO QCD corrections 0i and two-loop electroweak contributions i , is 

X{xt) = 1.481 ± 0.005th ± O.OOSexp, (2.4) 


where the hrst error comes from the remaining renormalisation scale and scheme 
uncertainties, as well as the theoretical error on the MS parameters due to the 
matching at the electroweak scale, while the second one corresponds to the combined 
experimental error on the top and W masses entering the ratio xt, and on the strong 
coupling as{Mz)- The central value and errors in (2.4) have been obtained using 


the MS couplings with full NNLO precision 34 - 3-loop running in the SM and 


2-loop matching at the weak scale (plus 4-loop QCD running of ag and 3-loop QCD 
matching in ag and yt) - and varying the renormalisation scale between Mt/2 and 
2Mt. The NLO EW correction has been included, using the result presented in [^, 
in order to eliminate the large EW renormalisation scheme dependence of the pure 
QCD result. See Appendix [A| for details about the different contributions to X{xt). 

The parameter Pc{X) summarises the charm contribution and is defined through 


Pc{X) = PS°(A) + 5P, 


6Pcu = 0.04 ±0.02, 


(2.5) 


with the long-distance contributions 5Pc^u calculated in 
lations could reduce the present error in this part 
given by 


10 


Future lattice calcu- 
35 . The short-distance part is 


pFix) = 4 


+ g^NNL 


where the functions result from QCD NLO js 


36 


( 2 . 6 ) 


and NNLO calculations 


. They also include complete two-loop electroweak contributions . The index 
distinguishes between the charged lepton flavours in the box diagrams. This 
distinction is irrelevant in the top contribution due to mt ^ mg but is relevant 
in the charm contribution as mr > me- The inclusion of NLO and NNLO QCD 
corrections have reduced considerably the large dependence on the renormalisation 
scale He (with Hc = 0{mc)) present in the leading order expressions for the charm 
contribution. The two-loop electroweak corrections on the other hand reduced the 
dependence on the dehnition of electroweak parameters. An excellent approximation 
for Pf^{X), including all these corrections, as a function of as{Mz) and rndme) is 
given in (50) of (see Appendix [B| . Using this formula for the most recent input 
parameters |37[|38 


A = 0.2252(9), mdme) = 1.279(13) GeV, ag{Mz) = 0.1185(6) (2.7) 


we find 

P^^{X) = 0.365 ±0.012. 


Adding the long distance contribution in (2.5) we finally find 


( 2 . 8 ) 


PdX) = 0.404 ± 0.024, 


(2.9) 
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where we have added the errors in quadratures. We will use this value in our 


numerical analysis. In obtaining the error in (2.9) we kept A fixed at its central 
value, as its error is very small and the strong dependence on A in P^^{X) 


IS 


canceled by other factors in the formula for the branching ratio as discussed above. 


2.2 Kj 




The branching ratio for Xl —)• 'k^vv in the SM is fully dominated by the diagrams 
with internal top exchanges, with the charm contribution well below 1%. It can be 
written then as follows 


39,40 


B{Kl —)• TT^Z/I/) = Kl • 


/ImAt 

V A5 


X{xt) 


where [TT 


kl = (2.231 ±0.013) • 10 


-10 


A 

0.225 


1 8 


( 2 . 10 ) 


( 2 . 11 ) 


We have summed over three neutrino flavours. An explicit derivation of (2.10) can 
be found in (^. Due to the absence of Pc{X) in (2.10), the theoretical uncertainties 
in B{Kl are due only to X{xt) and amount to about 1% at the level of 

the branching ratio. The main uncertainty then comes from ImAt, which is by far 
dominant with respect to the other parametric uncertainties due to kl and mt , with 
the latter present in X{xt). 


2.3 Experimental prospects 

Experimentally we have 

B{K+ ^ 7r+i/i7)e.p = (17.3l“;5) • 10-11, (2.12) 

and the 90% C.L. upper bound 

B{Kl vr°i/z^)exp < 2.6 • 10"® . (2.13) 


The prospects for improved measurements of B{K'^ —)• are very good. 

One should stress that already a measurement of this branching ratio with an accu¬ 
racy of 10% will give us a very important insight into the physics at short distance 
scales. Indeed the NA62 experiment at CERN 20,21 is aiming to reach this pre¬ 
cision, and it is expected to accumulate 100 SM events with a good signal over 
background figure by 2018. In order to achieve a 5% measurement of the branching 
ratio, which will be the next goal of NA62, more time is needed. The planned new 
experiment at Eermilab (ORKA) could in principle reach the accuracy of 5% 43 
Concerning Kl —?■ the KOTO experiment at J-PARC aims in the first 

step in measuring B{Kl —>■ at SM sensitivity and should provide interesting 

results around 2020 on this branching ratio 15,22 . There are also plans to measure 


this decay at CERN and one should hope that Eermilab will contribute to these 


^Unfortunately the US P5 committee did not recommend moving ahead with ORKA and it appears 
that the precision on B{K'^ —>■ n'^vv) will depend in the coming ten years entirely on the progress made 
by NA62. 
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\^k\ 

2.228(11) X 10-3 

38 


Fk 

156.1(11) MeV 

44 


SrpKs 

0.682(19) 

45 


Bk 

0.750(15) 

44 

46 


AMk 

0.5292(9) X 10-2 ps-^ 

38 


Fb, 

190.5(42) MeV 

44 


AMd 

0.507(4) ps-i 

45 


Fb. 

227.7(45) MeV 

44 


AM, 

17.761(22) ps-3 

45 


Fbs \l Bbs 

266(18) MeV 

44 


7 

(73.2«:3)° 

47 


e 

1.268(63) 


44 


IK.I 

0.2252(9) 

45 






Ar,/r, 

0.138(12) 

45 


Vcc 

1.87(76) 


48 



1.519(5) ps 

45 


Vet 

0.496(47) 


49 


BBs 

1.512(7) ps 

45 


Vtt 

0.5765(65) 


50 


OisiMz) 

0.1185(6) 

38 


Vb 

0.55(1) 

50 

51 


rndmc) 

1.279(13) GeV 

37 






Mt 

173.34(82) GeV 

52 







Table 1: Values of theoretical and experimental quantities used as input parameters. 

efforts in the next decade. The combination of K~^ —)• 'k^vv and Kl —)• is 

particularly powerful in testing NP. Assuming that NA62 and KOTO will reach the 
expected precision and the branching ratios on these decays will be at least as high 
as the ones predicted in the SM, these two decays are expected to be the superstars 
of flavour physics after 2018. 


3 CKM inputs from tree-level observables 


3.1 Determination of the branching ratios 


As discussed in the introduction, the CKM matrix can be determined by the tree- 
level measurements \ Vub\-, \ Vcb\-, |14s|, and the angle 7 of the UT. Although this is in 
principle the optimal strategy, it is currently marred by disagreements between the 
exclusive and inclusive determinations of both \Vub\ and \Vcb\ ~ for recent reviews 
see 53-55 . We proceed to present the latest results of both determinations, as well 


as our weighted average, with which we will give the SM predictions in what we call 
strategy A. 

The most recent exclusive determinations from lattice QCD form factors are 

44l[5^ 


iKftlexci = (3.72 ± 0.14) X 10-3, |Wb|e,,ci = (39.36 ± 0.75) x lO'^. 
The inclusive values are given by [44|[^ 

iKblinci = (4.40 ± 0.25) X 10-3, \Vcb\inci = (42.21 ± 0.78) x 10-3. 


(3.1) 

(3.2) 


We take a weighted average and scale the errors based on the resulting (specifi¬ 
cally, we follow the method advocated in [^), which gives 

iKfelavg = (3.88 ± 0.29) X 10-3, = (40.7 ± 1.4) x 10-3. (3.3) 









































































3 CKM inputs from tree-level observables 



Figure 1: Error budgets for the branching ratio observables B{K^ —>■ vr+uz/) and B{Kl —>■ 
7 r°z/u). The remaining parameters, which each contribute an error of less than 1%, are 
grouped into the “other” category. 


For the CKM angle 7 we take the current world average of direct measurements 

7 = (73.2+6;3)°. (3.4) 


Using this, together with |Fus| = A already given in (2.7), we can determine the full 
CKM matrix. 

In particular, we can determine the quantities Xt and Ac, which enter the ex¬ 
pressions for the branching ratios in ( 2 . 1 ) and ( 2 . 10 ), as functions of these input 
parameters. These expressions are: 


ReAi ~ iKfellUcbl cos7(l - 2A^) + (|Kb|" - \Va\'^)X ( 1 - — 
ImAt ~ iKftllUcbl sin7. 


ReAc ~ -A ( 1 - Y 


(3.5) 

(3.6) 

(3.7) 


which, with respect to their leading order in A, are accurate up to O(A^) corrections. 
The (exact) numerical values for ReA* and ImAt obtained from our three different 
choices of Vui, and Vcb in (3.1)-(3.3) are given in Tablel^ 

These expressions can then be directly inserted into an d (| 2 . 10 [ ) in order to 
determine the two branching ratios. Using our averages from (3.3) together with 
(3.4) gives 


B{K+ TT+uv) = (8.4 ± 1.0) X 10"^\ 
13{Kl 7 r°nn) = (3.4 ± 0 . 6 ) x lO'^F 


(3.8) 

(3.9) 


In Figure we show the error budgets for these two observables, and see that 
the CKM uncertainties dominate. In particular in the case of K~^ —)• we 

observe large uncertainties due to \ Vcb\ and 7 , while in the case of Kl —t the 

uncertainty due to \ Vub\ dominates but the ones from |Ucb| and 7 are also large. The 
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Figure 2: Dependence of the branching ratio observables B{K^ —)• 'k^vV) (left) and 
B{Kl — )■ TT^pp) (right) on the CKM parameter inputs \Vcb\, |14b| 0 ^ 6 ^ 7- The 95% C.L. 
bands in Vub, Vcb and 7 are shown in green, blue, and red, respectively. 


remaining parameters, which each contribute an error of less than 1%, are grouped 
into the “other” category. 

For convenience we give the following parametric expressions for the branching 
ratios in terms of the CKM inputs: 


B{K^ iT+pp) = (8.39 ± 0.30) x 10 


-11 



2.8 

7 

40.7 X 10-3 


73.2° 


1 0.74 


(3.10) 
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B{Kl 



(3.36 ±0.05) X 10“^b 


[ \Vub\ ] 

2 

[ 1 

2 

sin( 7 ) 

3.88 X 10-3 


40.7 X 10-3 


sin(73.2°) 

(3.11) 


The parametric relation for B{Ki —)• is exact, while for B{K'^ —)■ tt+uu) 

it gives an excellent approximation: for the large ranges 37 < \Vcb\ x 10^ < 45 
and 60° < 7 < 80° it is accurate to 1% and 0.5%, respectively. In the case of 
B{K~^ —)• TT'^ui') we have absorbed \Vub\ into the non-parametric error due to the 
weak dependence on it. The exact dependence of both branching ratios on |I4b|, 
\Vcb\ and 7 is shown in Figure]^ 

In order to obtain the values of ek-, S^Ks^ ^^s,d and of the branching ratios 
for Bg^d ''^6 use the known expressions collected in 16 , together with 

the parameters listed in Table The “bar” on the Bg —)• fjT branching ratio, 
B{Bs —)■ denotes an average over the two mass-eigenstates, as measured by 

experiment, rather than an average over the two flavour-states, which differs in the 


Bg system 5^60 


In Table we show the results for the —>■ Tr~^iyu and iFn —>■ branching 

ratios and other observables, for three choices of the pair (|I4b|) Wcb\) correspond¬ 
ing to the exclusive determination (3.1), the inclusive determination (3.2) and our 
average (3.3). We use (|3.4[) for 7 in each case. We observe: 


• The uncertainty in B{K'^ —)• vr+nn) amounts to more than 10% and has to 
be decreased to compete with future NA62 measurements, but finding this 
branching ratio in the ballpark of 15 x 10“^^ would clearly indicate NP at 
work. 


• On the other hand, consistency with B{Bg —?■ would imply the —?• 

branching ratio to be in the ballpark of 7 x 10“^^. In such a case the 
search for NP in this decay will be a real challenge and the simultaneous 
measurement of —)• will be crucial. 


• The values of S^Ks typically above the data but only in the case of the 
inclusive determinations of both \Vcb\ and |I4ft| is a new CP phase required. 

• The accuracy on the SM prediction for AM^ and AM^ is far from being 
satisfactory. Yet, the prospects of improving the accuracy by a factor of two 
to three in this decade are good. 


3.2 Correlations between observables 

Correlations between and Bg 

From inspection of the formulae for the branching ratios for —)• and 

Bg^d l^~ 1 each of which in particular depends on \ Vcb\-, we derive the following 
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exclusive 

inclusive 

average 

measured 

BR(iC+ ^ n+uu) [10-11] 

7 fi9+0-69 
' •'J^-0.70 

q qq+0.89 
'9-0,92 

s qq+i-oe 

1 7 0+11.5 
' •'^-10.5 

BR{Kl ttVp) [10-11] 

0 qq+0.30 
^•^^-0.35 

4 64+®'®3 

q qc+0.60 

'J-'J'J-o.ei 

< 2600 

BR(i?s —t g~^g~) [10-®] 

q 1 0+0.18 

3 66+®'^i 
O.OO_0 20 

3.40l®;f7 

2.8 ±0.7 

BR(Rd ^ g+g-) [lO-i®] 

1 OO"*"®-!! 

1 17+0-14 
-‘-•■‘-'-0.14 

1.081®:J3 


[10-3] 

1 qR+0.25 
l.yU_o.27 

0 ^7/1+0.36 
'^-0.38 

Q qo+0.35 
^•^^-0.36 

2.228 ±0.011 

CSM 

‘^gKs 

o-74lS:SI 

0.807SS 

n 7c:+0.05 
' ^-0.05 

0.682 ±0.019 

AM, [ps-1] 

le.io^li 

18.64711 

17 44+2-74 

i 1 .0+_2,58 

17.761 ± 0.022 

AM, [ps-i] 

'9-J^-0.09 

0.60;®;}J 


0.510 ±0.003 

Im(At) [10-4] 

f 40+®-®'^ 

1 70+0.12 

I- 'O-0.13 

1 ^1+0.13 
-‘-•^■‘--0.14 

— 

Re{\t) [10-4] 

—2 QQ+®-i® 

—3 3q+®-24 

O.Oa_o 23 

—3 20+®'^® 

O.ZU-0 29 

— 

Rb 

n 41 +0.02 
'^•^■‘--0.02 

n 41^+0.03 
U-+^-0.03 

n 41 +0-03 
'^• + ■1-0.03 

— 


Table 2: Values of —>■ B{Kl —>■ vr^nn) and of other observables within the 

SM for the three choices of |14b| and \ Vcb\ following strategy A as discussed in the text. 


approximate relations 


i3(K+ ^ TT+nn) = (8.39 ± 0.58) x 10 


-11 


7 

73.2° J 


1 0.81 


B{Bg —>■ p+/i ) 

1.42 

■ 227.7 MeV 

3.4 X 10-® 


i Fb, \ 


2.84 


, (3.12) 


B{K+ TT+nn) = (8.41 ± 0.77) x 10” 


11 

B{Bg — >■ g~^g ) 

0.74 

■227.7MeV 


3.4 X 10-® 


i Fb. \ 


1.48 


B{Bd — )• p+p ) 

0.72 

'190.5 MeV 

1.08 X 10-1® 


i Fb, \ 


1.44 


(3.13) 


Note that both relations are independent of |14fe| and (3.12) depends only on 7 . In 
particular the correlation ( 3.12[ ) should be of interest in the coming years due to the 
measurement of K~^ —)• by NA62, of Bg —)• gAg by LHCb and CMS and of 7 

by LHCb. Moreover the last factor should also be improved by lattice QCD. 

In the left panel of Figure]^ we show the correlation between —?■ 7r'’'nn and 

Bg —)■ p+/i for different fixed values of 7 . The dashed regions correspond to a 68% 
C.L. that results from including the uncertainties on all the other input parameters, 
whereas the inner hlled regions are a result of only including the uncertainties of 


\Vub\, \Vcb\ (we use the averages in (3.3)), and |14s|. 

It should be noticed that the present experimental determination of B{Bg —?• 
g~) is slightly lower than the SM prediction, and the agreement between the 
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Figure 3: Left panel: correlation of B{Bs —)■ in~) versus B{K^ —)■ tt+i/i?) for fixed 
values 0 / 7 . Right panel: correlation of B{Kl —)• vr^i/p) versus B{K^ —)■ vr+z/P) for fixed 
values of fd. In both plots the dashed regions correspond to a 68% C.L. resulting from the 
uncertainties on all other inputs, while the inner filled regions result from including only 
the uncertainties from the remaining CKM inputs of strategy A. 


SM and the data can be improved by lowering \Vcb\ to values in the ballpark of 
its present exclusive determinations. But in this case, as can be seen already from 
Table the SM predictions for both —)• vr+nP) and ek are also reduced. It 

can be useful to express B{K~^ —)• vr+nP) as a function of ek, in a way similar to 
( |3T2l ) and ( [3T3l ), in order to make the correlation between them explicit. One has 


B{K"^ vr+nP) = (8.39 ± 1.11) x 10 


-11 


\£k\ 


2.23 X 10 


-3 


1.07 


r 7 

- 0.11 

[ \Vub\ 1 

L73.2°J 


3.88 X 10-3 


1 -0.95 


(3.14) 


We do not write explicitly the dependence on the hadronic quantities, since here 
more parameters are involved. The uncertainty here comes mainly from g^c and get, 
while the ones due to Fk are smaller than the corresponding ones in the Bg^d meson 
systems. It is evident from this formula that a reduction of B{K~^ —)• 7r'’'nP) implies 
also a reduction of ek- 


The correlations in (3.12), (3.13) and (3.14) result from the fact that it is possible, 
by taking suitable powers of the B^ d branching ratios, to eliminate the 

dependence on \Vcb\-, while the one-loop functions X, Y, and S are hxed by the 
top mass in the SM. Both correlations could be broken already in models with 
constrained MFV (CMFV) in which the modiheations of the functions X and Y are 
generally different. In general MFV models new scalar operators could additionally 
contribute to Bg^d 6%, modifying also the factors involving the weak decay 
constants. Therefore, these correlations are strictly valid only in the SM, and their 
violation would not necessarily rule out (C)MFV. 
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7T~^uu and Kl in Minimal Flavonr Violation 

In models of NP with Minimal Flavour Violation (MFV) there are no flavour¬ 


changing interactions beyond those generated from the SM Yukawa couplings 61 
For K~^ —)• TT+nn and Ki ^ this restricts the operators that can contribute in 
most NP models to just the operator already dominant in the SM, (sd)y_^(nn)y_^, 
making MFV equivalent to Constrained MFV (CMFV) in this case. Therefore 
in MFV the value of X{xt) can be shifted but must stay real. Defining for conve¬ 
nience 


B. = 


B{K^ —7- 

k+(1-|-Aem) 


Bl = 


B{Kl —)• 


kl 


(3.15) 


with K+ and kl given in 
the correlation 

B+ = Bl + 


and (2.11), respectively, we have in the case of MFV 


“ y ) sgn(V(xO)Fc(V) 


(3.16) 


which was first given in 63,64 . Note that in the SM the sign of X{xt) is positive. 
Recalling the relation 

^=.-cot,3(l-M , (3.17) 


accurate to terms, and solving for /3 then gives 


cot j3 = 


1 




B+-Bl 


Bl 


-'-Y 


sgn(V(xt))Pc(Y) 


(3.18) 


In deriving (3.18) from ( |3.16 ) and (3.17) one finds that for X{xt) > 0 this 
solution is unique, while for X(xt) < 0 a second solution with a minus sign in front 
of first square root is allowed [^. However this solution is excluded if we require 
both branching ratios to be larger than 10 “^^ and we will not consider it here. 

In the SM and CMFV we have to a very good approximation the relation |63|64| 

SipKs = sin 2/3, 


(3.19) 


which is only spoiled by possible penguin enhancements in the B^ —)• J/'ipKg 


mode 65 
K+ 


Thus (3.19) together with (3.18) give a triple correlation between 
and S^Ks fbe SM and CMFV. 


Kl 

As demonstrated in the earlier parts of this section, the branching ratios for 
K+ —)■ 7 r“*'nn and Kl —?• still contain significant parametric uncertainties due 

to the uncertainties in \ Vcb\-, \Vub\ and 7 , and to a lesser extent in mt- It is therefore 
remarkable that within the SM all these uncertainties practically cancel out in this 


triple correlation 63 . Moreover, this property turns out to be true for all models 


with constrained MFV 64 


We note that the main nncertainty in (3.18) resides in Pc, as the uncertainty in 


A is very small. We stress that this relation is practically immune to any variation of 
the function X{xt) within MFV models. This means that once B{K^ —)• and 

SipKs wifi be precisely measured we will know the unique value of B{Kl —>• 
within CMFV models. This relation is analogous to the one between B{Bs^d 
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MFV relation — present constraints 


" IKislimci/llrtliiici central value 

• |l'ti6|exci/|tc6|exci Central value 

□ 95% CL using |Vb,|.vs/|Kf,Uvg 
95% CL using S^-ks ^ sin 2/3 
68% CL A''*" —^ measurement 


MFV relation — future constraint£^^2025^ 



BR(A+ 


Figure 4: The MFV relation between —)• vr+uz/ and Kl —)■ vr^z/u using S^Ks — sin 2/3 

versus using the various tree-level inputs of | Vfe/Kifel and 7 (see text). In the left panel we 
show situation from current constraints, and in the right panel the possible situation in 
the following decade, including 10% precision on the two branching ratios, for illustration. 


and AMg^d 66 , where the present knowledge of AMg^d together with the 
future precise value of B[Bg —)• pV) will allow us to uniquely predict the branching 
ratio B{Bd —)• pVpi~) in CMFV models well ahead of its precise direct measurement. 

In the right panel of Figure]^ we show the correlation between —)• -n^vU and 

Kl ^ TT^uu for different fixed values of /3 {S.^Ks)- Tbe dashed regions correspond 
to a 68 % C.L. that results from including the uncertainties on all the other input 
parameters, whereas the inner filled regions are a result of only including the un¬ 
certainties of \Vcb\ (we use the average in (3.3)), 7 (as given in (3.4)) and |V, 


m 


(2.7). We observe that in the latter case the dependence on the remaining CKM 


parameters, for fixed (3, is indeed minimal. 


It is also possible to express the ratio in (3.17) as 


ReAt 

ImA^ 


I- ^ 
^ 2 


sm7 


cos 7 



Fc6 


Vuh 


(3.20) 


i.e. in terms of the tree-level CKM inputs discussed in this section, which are gen¬ 
erally assumed to be free of NP effects. We note that in MFV also S.,j,Ks is not 
affected by NP and is more accurately determined than 7 . On the other hand, 
there is a class of models - e.g. models with a U{2)^ flavour symmetry 67 - where 


the correlation with S'l/iXs is no longer true, while (3.16) and the generic relation 
(3.20) still hold. 


While the virtue of the correlation (3.18) is its very weak dependence on the 


CKM parameters, the correlation (3.16) together with (3.20) shares partly this 


property as it depends only on the ratio \Vcb/Vub\, equivalent to Rb, and not on 
\Vub\ and I Feb I separately. As we can see from the values of Rb given in Table 
this avoids some of the trouble with exclusive versus inclusive determinations, as 
the ratio of purely exclusive or inclusive determinations, as well as their weighted 
average, results in less variation - i.e. only 5% among the cases considered. Note 
that combining exclusive \Vub\ with inclusive |F;b|) for example, gives a greater 
variation. 
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In the left panel of Figure we compare the MFV relation for various values 
of iVcft/Fufel, including ala C.L. region corresponding to our weighted averages. 
We also include for comparison the relation corresponding to the current S.^Ks 
measurement, which is seen to be more accurate. In the right panel we repeat this 
comparison for possible future constraints in the following decade, where our choice 
of errors are based on those collected in . For the branching ratios of —)• 7r'*'nn 

and Kl —>• we assume a 10% precision relative to the SM predictions. Though 

this is the realistic target set by the NA62 experiment for K'^ —)• 7r'’'nn, the KOTO 
experiment will likely not reach such a precision for Kl —)• We observe that 

in this possible sketch of the future, the two decays under consideration have the 
potential to probe MFV and/or a U{2)^ symmetry. 


4 CKM inputs from loop-level observables 

A different approach is to assume that there are no relevant NP contributions to 


all the quantities listed in (1.2), so that we can use them together with the precise 
value of 1 14 s I to determine the best values of 


/3, 


|%cb|, \Vub\, \Vtd\, 


IK 


ts I ? 


(4.1) 


and predict the branching ratios for —)■ Kl —)• and Bg^d 


Clearly the absence of NP effects in all the loop observables (1.2) requires the SM 


to be valid up to a reasonably high energy scale, which is a stronger assumption 
with respect to the one of strategy A, where only tree-level determinations of CKM 
parameters were assumed to be free of NP effects. We call this approach strategy B. 


The relevant SM expressions can be found in 16 and in particular in |30|, where 


this strategy has been used to determine the correlation between the values of | 14 b| 
and | 14 b| with the non-perturbative parameters relevant for AMg^d- As the precision 
on these parameters resulting from QCD lattice calculations is improving, and the 
value of Bki relevant for ek, has been known precisely already for some timej^we 


can now use these formulae to extract the values listed in (4.1). 


At least three independent observables among the four listed in (1.2) have to 


be used in order to fix the three free parameters of the CKM matrix besides |14s|- 
For illustration, we present here the strategy which allows us to determine the 


parameters in (4.1) with high precision with the minimal number of measurements. 


Schematically this procedure can be described in two steps: 

• Step 1: The unitarity triangle can be determined from the experimental values 
of S^Ks = sin 2/3 and the mass ratio 


AMd niBa 1 iKdl" 


AMg 




Fb, \ Bb, 


(4.2) 


FBayBBd 

Using the following very accurate expressions for \Vtd\ and |Ks|) 


\Vtd\ ^ MVcblRt, 


IK 


ts 


l + y(l-24?tcos/3) ) iKfel, 


(4.3) 


®We use Bji = 0.750(15) which takes into account the values obtained by lattice QCD 44 and large 


N approach 46 . 
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where Rt is one of the sides of the UT, and solving (4.2) for Rt one gets 

I 


Rt ^ 


_ I ms, 

\Vus\\l AMs rUBa 


V AMs ruB, 2 


(4.4) 


where the dots indicate terms of order 0 ( 114 s| 14 spAMrf/AMs). 

With one side, Rt, and one angle, (3, known, the full unitarity triangle is 
determined by means of purely geometrical relations. In particular one has 


Rb 


1 - A72 

A 


Vcb 




2i?t cos/3, 


cot 7 


1 — Rt cos /3 
Rt sin /3 ’ 


(4.5) 


and the apex {g,r]) of the triangle is given by 


Q = 1 — Rt cos /3, T] = Rt sin /3 . 


(4.6) 


The very precise value of Rt obtained through AM^/AMs therefore allows a 
very precise determination of 7 . 

It should be emphasised that the UT constructed in this manner is universal 
for all CMFV models as the box function S does not enter the expressions 
used in these two steps [^. Moreover this determination is independent of 

• Step 2 : The measured value of {ekI then allows us to determine the optimal 
value of \Vcb\- Indeed we have 




= K, 


M^ 

' 6^/27r2AMi^ 


BK\V,b\^X^Rt sin/3, 


X {\Vcb\‘^Rt cos f3r]ttSo{xt) + T]ctSo{xc, xt) - ticcXc) 


(4.7) 


where xt = Sq is the well known SM box function as defined e.g. 

and Ks = 0.94 ± 0.02 


m 


16 


68 , 6 ^. With Rt known from (4.3) and /3 


determined from S^Ks^ fbe only unknown in (4.7) is |I4b|- Having found |I4b|, 
Rt,, and Rt, also |I4f,|, \Vtd\ and \ Vts\ are determined by the previous relations 
in a straightforward way. 

Alternatively we can also determine |I4b| by using separately AM* oc |V)sp and 
AMrf oc \Vtd\^ instead of ek, as both are proportional to |I 4 fep via the expressions 
given in ( |4.4| ). In principle it is also possible to determine the CKM matrix from 
AMd, AMs, and ek, but the precision in this case will be rather limited, due to the 
absence of the strong constraint on /3 from S,pKs ■ The best accuracy is obtained by 
performing a simultaneous fit to all the four observables AMd, AMs, S^pj^s ek- 
In Table we give the results of fits for the CKM matrix elements using dif¬ 
ferent combinations of the inputs discussed in the steps above. The values for the 
experimental and lattice observables used as inputs are listed in Table The fits 
were performed using a Bayesian statistical approach: uncorrelated Gaussian priors 
were chosen for each of the input parameters and the posterior distribution was 
sampled using Markov Chain Monte Carlo with the help of the Bayesian Analysis 
Toolkit 70 . A direct minimisation of the yielding identical results, has also 
been performed as a check. We observe that using \£k\ in step 2 gives a more 
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{{Ek], AMs, S'l/.XsfsM 

{AMd, AMs, +ifg}sM 

{\ek\, AMd, AMs, +ifs}sM 

|Uc.| [10-3] 

42.594}:" 

41 . 304 I 47 

42.354}:}3 

\Vub\ [10-3] 

3.6240:}3 

q Ki +0.27 
'^•^-‘--0.25 

'^•^-^-0.14 

m [10-3] 

8.9640:1 

8.684g:^^ 

8.954^:!! 

IU.I [10-3] 

41.794}:^ 

40.5241^^ 

41.554}:?! 

B{K^ —)■ Mvv) [10-3^] 

9.184^:}? 

8.394}:!? 

9-084g:!^ 

B{Kl n°iyi>) [10-“] 

0 n 1+0.33 
'^•'^-‘--0.29 

2.6642:^! 

2 Qg+0'32 

^.yo_o.28 

B{B, ^ Mk) [10-®] 

0 f;Q+0.30 

O.uy_o.26 

3.4642:!^ 

0 (7 1+0.27 
'^•^^-0.24 

B{Bd ^ Mm) [10-13] 

1 nn+O-OS 
-‘-•'^^-0.08 

1 02+^-^’^ 
-*^•^^-0.15 

-‘-•^^-0.08 

Ini(Af) [10-^] 

-*-•^'^-0.07 

1 qc;-l-0-20 

-*-•^^-0.07 

Re(Ai) [10-^] 


q 9C+0.40 
'^•^'-’-0.45 

'^•^'^-0.18 


Table 3: Results of the fit to the CKM matrix elements for various combinations of inputs 
as detailed in strategy B, and the corresponding observable predictions. 


precise result for \Vcb\ than the alternative of using and AM^ separately, as 

well as favouring a higher central value. The most accurate determination (given in 
the last column of the table), follows from including all inputs. The corresponding 
CKM matrix elements of interest are: 

\Vub\ = (3.61 ± 0.14) X 10-3, \Vcb\ = (42.4 ± 1.2) x lO'^, 

\Vtd\ = (8.94 ± 0.27) X 10-3, \Vts\ = (41.6 ± 1.2) x 10-3. (4.8) 

For completeness, we give here the sides of the UT as determined from onr full 
fit, that read 


Rb = 0.368 ±0.013, 


= 0.937 ± 0.032, 
while its angles are 

a = (89.0 ±5.0)°, /3 = (21.5 ±0.8)° 

and its apex 

^ = 0.129 ±0.030, r/= 0.344 ± 0.017. 


7 = (69.5 ±5.0)°, 


(4.9) 

(4.10) 

(4.11) 


The precision on Rt, 7 and |14b| using the above strategy is already impressive, 
and will continue to improve with new lattice results. Using for instance the im¬ 
proved error estimates for ^ and fs^ \lB bs from 31 (keeping the central values 
from [44]) we find the very precise results: 


lUcbl = (42.0 ±0.9) X 10-3, 7 = (70.8 ±2.3)°, = 0.945 ± 0.015. (4.12) 

In Fignrej^we show the htted ranges for \ Vub\ and |14fe| and compare them with 
the inclusive, exclusive and our averaged values in (3.1)-(3.3). We distinguish be¬ 
tween three different cases: the bine area corresponds to the fitted range of \ Vub\ and 
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0.044 


0.042 


^ 0.040 


0.038 


0.036 


0.0030 0.0032 0.0034 0.0036 0.0038 0.0040 0.0042 0.0044 0.0046 0.0048 

\Vub\ 

Figure 5: Comparison of 68% C.L. regions for \Vuh\ and \Vcb\ in strategy B for various 
combinations of inputs versus their reported inclusive and exclusive values, and our aver¬ 
ages of these, as considered in strategy A. 

I Feb I determined by \£k\-, AMc^/AM^ and S^Ks'i for the green area AM^, AM* and 
S^Ks r^re used as inputs and the red area combines both and uses \£k\, AM^^, AMg 
and SijjKs s-s inputs. As noted earlier, one can see that especially le^"! favours large 
values of |Fcb|, around the inclusive value, while the rather small \Vub\-, around the 
exclusive value, is favoured by S^Ks- 

It is interesting to compare these results with the indirect fits performed by 
and CKMfitter [^, which give 

UTfit: \Vub\ = (3.63 ± 0.12) x 10■^ \Vcb\ = (41.7 ± 0.56) x 10“^ , (4.13) 

CKMfitter: \Vub\ = (3.551°;^^) x lO’^, |F,b| = (41.17;°;?°) x 10"° . (4.14) 

They are in very good agreement with our results. We note however, that these two 
groups included in their analyses the information from tree level decays, which we 
have decided not to include in our strategy B because of the discrepancies between 
inclusive and exclusive determinations of |Fife| and \Vcb\- Moreover, we also did not 
use the tree-level determination of 7 contrary to these two groups. 

Having determined the full CKM matrix in this manner, predictions for rare 
decays branching ratios can be made. These are collected in the last four rows in 
Table and again the most precise are the ones in the last column so that our final 
results for the four branching ratios are: 



= (9.11 ±0.72) X 10"^; 

(4.15) 

B{Kl — )• 

= (3.00 ±0.31) X 10"^; 

(4.16) 

B{Bg — )• ) 

= (3.66 ±0.26) X 10“°, 

(4.17) 

B{B(i —)• pA g ) 

= (1.09 ±0.08) X 10“^°. 

(4.18) 


UTfit 27 


inclusive measurements 


«.-- 

I our average 


exclusive measurements 

Strategy B: assuming SM in loops 

I-1 i\€K\,AMJAM„ S^kAsm @ 68% CL 

I-1 {AMj, AM,, Sv.A'JsM @ 68% CL 

I-1 {kA-l. AM,. AM,, S,pK,}sM @ 68% CL 
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Figure 6: Comparison of 68% C.L. regions for B{Kl —)• 7r°i/z/) and B{Bs —)■ pCja ) versus 
B{K^ —)■ 'k^uD), using different inputs from both strategy A and B to fix the CKM matrix 


In (4.12) we used the new lattice error estimates from 31 for a “sneak preview” of 


how the CKM fit in strategy B will improve once the full results will be available. 
Using these results for the observable predictions listed above will likewise lead to 
reduced uncertainties: SB{K~^ —>■ tt+uu) = 0.65, 6B{Kl —)• = 0.28, 6B{Bs —)> 

pL~^ia~) = 0.22 and SB{Bd —)■ pA) = 0.07. 

As a comparison, using instead the fit results of (4.13[) and (4.14), one gets 


UTfit: B(A+ ^ tt^uP) 

= {SMflit) X 10-11, 

(4.19) 

B{Kl —)• AvP) 

= (2.93 ±0.25) X 10-11, 

(4.20) 

CKMfitter: B{K^ —)• tt^vP) 

= (8.17;0;°J) X 10-11, 

(4.21) 

B{Kl —)• AvP) 

= (2.65;°:i) X 10-11. 

(4.22) 


It is also interesting to compare the results in ( |4.17 ), ( |4.18 ) with the most recent 
prediction in the SM 71 , with which our SM results are in perfect agreement]^ and 
with the most recent averages from the combined analysis of CMS and LHCb 
that read 


B{B, ^ p+p-) = (2.8l°:^) X 10“®, (4.23) 

B{Bd ^ p^p-) = (3.9li:®) X 10-1°. (4.24) 

Note that the SM value of B{Bs —)■ p~^p~) is outside one sigma range of the exper¬ 
imental value. 

In Figure the correlations of B{Kl —)■ AuP) and B{Bs —)■ p~^p~) versus 
B{K~^ —)• TT^uP) are shown, comparing the best result of strategy B, which in¬ 
cludes all of the available inputs, with the inclusive, exclusive and average cases 
of strategy A. We observe that the inclusive case of strategy A is very similar to 
strategy B for 77+ —)■ -/r+nn and Bg —)■ p'^p~, as both have little sensitivity to |I4b|, 

®This is not surprising as these authors used the inclusive determination of \Vcb\ that is very close to 
the value determined by us. 














































5 The ratio s' je in the Standard Model 


20 


whereas Ki^ —?• which has a stronger \ Vub\ dependence, can differentiate them. 

In both plots our average for \Vub\ and \ Vcb\ is seen to also pick the middle ground 

for these observables. _ 

Evidently, the present experimental value for B{Bs —)• in (4.23) would 

favour the exclusive determination of \Vcb\ and a value of B{K~^ —)• vr+nn) in the 
ballpark of 7 x 10“^^ rather than 9 x 10“^^. But then also the value of \sk\ would 
be below the data. It appears then that unless the experimental value for B{Bs —)• 
moves up by 20% in the coming years, the SM will face some tensions in this 
sector of flavour physics. 

It is instructive to recall the following formula [3 ,73 
pendence of B{K~^ —)■ vr+z^T') on Rt, (3 and Vcb'- 


that summarises the de- 


B{K^ 


iVcbl'^Xixtf 


+ 


1-^ 

2 


sin2 /5 


(1 - A2/2)2 


Rt cos (3 + 


\Vcb\^Xixt)) 


(4.25) 


This can be considered as the fundamental formula for a correlation between B{K~^ —)• 
and any observable used to determine i?t, and is valid also in all models 




with MFV where X{xt) is replaced by a real function X. When this formula was pro¬ 
posed, it contained significant uncertainties in Rt determined through AM^/AMg, 
in Pc{X) known only at NLO, in and in \ Vcb\- The first three uncertainties have 
been significantly reduced since then. Moreover, the improved knowledge of the 
non-perturbative parameters entering ek and AMg^d allows now within the SM to 
determine \ Vcb\ rather precisely. We stress that in other models with MFV the latter 
determination will depend on the NP contributions to ek and AMg^d which modify 
the function S. An analysis of this issue is presented in 
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Finally when 7 from tree-level decays will be precisely measured, Rt will be 
determined solely by (3 and 7 , 


Rt = 


sm7 


sm{(3 -I- 7 )' 


(4.26) 


and the dependence of B{K^ —?■ on 7 can be directly read off (4.25). 


5 The ratio e'/e in the Standard Model 

The ratio e'/e measures the size of direct CP violation in —)• vr-Tr relative to 

the indirect CP violation described by ek- In the SM e' is governed by QCD 
penguins, but receives also an important destructively interfering contribution from 
electroweak penguins that is generally much more sensitive to NP than the QCD 
contribution. 

The ratio e'/e is measured to be 

Re(e7e) = (16.5 ± 2.6) x 10"^, (5.1) 

and the imaginary part of e'/e is negligible so that we will just write e'/e in all 
formulae below. 
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This result constitutes in principle a strong constraint on theory. However, the 
difficulty in making predictions for e'je within the SM and its extensions is the 
strong cancellation of QCD penguin contributions and electroweak penguin contri¬ 
butions to this ratio. In the SM QCD penguins give a positive contribution, while 
the electroweak penguins a negative one. In order to obtain a useful prediction for 
s'/e in the SM the corresponding hadronic parameters Bg ' and Bg ' ' have to be 
known precisely. 

In the large N limit one has = 1 

1/N corrections in 80 


77 


79 


While the study of 


( 3 / 2 ) 

indicated that Hg ' < 1, no conclusive result has been 

obtained for Fortunately, very recently significant progress on both 

( 3 / 2 ) 

and Bg ' ^ has been made by lattice QCD simulations and in the context of the large 


N approach. Indeed, from the results of the RBC-UKQCD collaboration |81, 82 
one can extract 


B. 


( 1 / 2 ) 


= 0.57 ±0.19, 


(3/2) 


(lattice QCD) (5.2) 


83 


as shown in Appendix in the case of Hg , and in 
other hand, the very recent analysis in the large-approac’ 
a conservative upper bound on both and which reads 


84 


for B, 


( 1 / 2 ) 

6 


On the 


< 1 . 


1 in |84 allows to derive 


(5.3) 


(large-A^) 


(o i‘2'\ 

Moreover, one finds Bg ' {me) = 0.80 ± 0.10 in good agreement with (5.2). The 


result for B, 


( 1 / 2 ) ■ 


g is less precise but there is a strong indication that Bg 


(l/2l 


< B 


(3/2) 


(3/2) 


also in agreement with (5.2). We refer to 84 for further arguments why is 


expected to be smaller than Bg 

The most recent analysis oie' je has been given in 83 . Using the results in (5.2) 
and determining the remaining contributions to e' je by imposing the agreement of 
the SM with CP-conserving data one find^ 83 

(5.4) 


Re(e7e) = (1.9 ±4.5) x 10 


-4 


significantly below the experimental value in ( |5.1[ ). This result differs by roughly 
3(7 from the data, but, as stressed in [^, larger values can be obtained if only the 
absolute large N upper bound on both parameters in (5.3) is used. Yet, as found 


there and confirmed here by us, even with more generous values of ImAj the SM 
has serious difficulty in describing the data for e'/e. 

In spite of this it is of interest to study the correlation of e' je with —)• tt^vv 

in t he SM as this correlation has been already studied in various extensions of the 

SM this correlation depends only on the values of 
and the CKM parameters which we determined in the previous sections 


SM 
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and Bo 


85 

75T 


using strategies A and B. 

All the relevant details on e' je within the SM including the relevant references 
Here we collect only the relevant information necessary to perform 


are given m 


83 


the numerical analysis. The basic analytic formula for e' je reads 83 


= ImAt • Fe'{xt), 


(5.5) 


SM 


^To this end Im A* = (1.4 ± 0.1) x 10 ^ has been used. 
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where 


F^t{xt) = Po + Px Xo{xt) + Py Yo{xt) + Pz Zo{xt) + Pe Eo{xt) , (5.6) 


with the first term dominated by QCD-penguin contributions, the next three terms 
by electroweak penguin contributions, and the last term being totally negligible. 
The xt dependent functions have been collected in Appendix [C} 

The coefficients Pi are given in terms of the non-perturbative parameters Rq and 
Rs defined in (5.8) as follows: 


Pi = rP + rP Re + rP Pg • (5.7) 

The coefficients rP, rP and rP comprise information on the Wilson-coefficient 
functions of the AS" = 1 weak effective Hamiltonian at NLO. Their numerical values 
are given in the NDR renormalisation scheme for /x = nic and three values of as{Mz) 
in Table in Appendix 

The parameters Rq and Pg are directly related to the parameters and 

( 3 / 2 ) 

Pg ' ' representing the hadronic matrix elements of Qq and Qg, respectively. They 
are defined as 


Pe = P( 


( 1 / 2 ) 

6 


Pg = B, 


(3/2) 

8 


114.54 MeV 1^ 
msirric) + mdirric) \ 

114.54 MeV 
ms{mc) + mdirric) _ 


(5.8) 

(5.9) 


(1/2) (3/2) 

We stress that both PI ' ^ and Pg ' ' depend very weakly on the renormalisation 


scale 91 


In Figure]^ we show the correlation between e^/e and Ke —)• in the SM. 

The central value from the RBC-UKQCD collaboration in (5.2) has been used for 

(3/2) (1/2) 

Pg h The different colours correspond to different choices of the parameter Pg h 



(blue), 

(5.10) 


(green), 

(5.11) 


(red). 

(5.12) 


The first choice is motivated by the upper limit from large N approach in (5.3), 
although the bound Bq ' < Pg Ms violated, and gives an idea of the largest 
possible values of e' je attainable in the SM. The second choice assumes that B^"^^ = 
Pg Ms saturating the previous bound. Finally, the third choice uses the central 
values (5.2) from the RBC-UKQCD collaboration for both Pg^^^^ and bP‘^\ We 
observe that even for the first choice of Pg and Pg ' the ratio e'/e in the SM 
is below the data, and only for the largest values of ImA^ it is within 2a from the 
central experimental value. For such values also the branching ratio for Kl 


TT^VV 


is largest. 
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Figure 7: Correlation ofB{KL —)■ 7r°z/P) versus \e'/e\ for fixed values of = 0.57 (red), 
0.76 (green), and 1.00 (blue). The hatched regions correspond to a 68% C.L. resulting 
from the uncertainties on all other inputs, for strategy A using our average values of 114x61 
and 114:61 7 and strategy B. The yellow band shows the experimental result at la. 


: = 0.57 


0.76 



r//i Strategy A - 68% C.L. 
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6 Summary and outlook 


In this paper we have performed a new analysis of the rare decays —)• and 

Kl —)• within the SM. The prime motivations for this study were: 

• The start of the NA62 experiment that should in the coming years reach a 
precision of 10% relative to the SM prediction for B{K^ —)> -k^vu). 

• The soon to improve value of for which preliminary error estimates are 
already given in [^, which will allow a much more precise determination of 
the elements of the CKM matrix, in particular of the angle 7 , 114x61 and |146|) 
without the use of present tree-level determinations of these parameters that 
are presently subject to significant uncertainties. 

• The observation of the correlation between B{K'^ —?■ B{Bs —)• pAix~) 

and 7 within the SM that only weakly depends on |Fc 6 | • This correlation should 
be of interest in particular for CERN experimentalists who in the coming years 
will significantly improve the measurements on these three quantities. 


Our main results are illustrated with several plots in Sections and Our analy¬ 
sis demonstrates that in the coming years the SM will undergo an unprecedented test 
due to the measurements of the rates for the decays —)• and Bg —)• p~ 

and improved determinations of the CKM parameters either through the strategies 
A or B, accompanied by improved lattice QCD calculations of the relevant non- 
perturbative parameters. Around 2020 these studies will be enriched through precise 

K{K* 


measurements of the rates for Kl 


Aw, 


Bd-f p'^p and Bd 

41/2) 




(n /O'^ 

Also improved knowledge of the parameters and B^ , accompanied 

with improved values of CKM parameters, will allow a more precise prediction for 
the important ratio e'/e. Calculating this ratio using strategies A and B we find. 
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in accordance with the recent analysis in 83 , that the SM prediction for e'/e is 
significantly below the data leaving large room for NP contributions. A recent 
analysis of e' je in simplified NP models has shown which NP models could move 
the theory prediction for e'/e to agree with data 


89 


Needless to say, in order to 
be sure that the SM indeed fails in the description of data, a big effort in clarifying 
various uncertainties will be required, as discussed in . 

It should be observed that the agreement of the SM prediction for B{Bs —t- 
with the data can be significantly improved by lowering |I4fe| to the values in 
the ballpark of its present exclusive determinations using lattice QCD form factors. 
But then automatically ek is found significantly below the data. Interestingly in 
this case B{K~^ —)■ is also predicted to be in the ballpark of 7 x 10“^^, that 

is more than a factor of two below its present experimental average. No doubt, the 
coming years will be exceptional for quark flavour physics. 
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A Expression for Xt 

The loop function Xt of ( |2.4[ ) can be written as 

/y Qj 

X{xt) = Xo{xt) + -^Xi{xt) + —XEw{xt), (A.l) 

47r 47r 

where Xq is the leading order result, and Xi, Vew are the NLO QCD and EW 
corrections, respectively. The coupling constants Og and a, as well as the parameter 
Xt = nil/m^ = 2yllg\^ have to be evaluated at a given renormalisation scale 
ii^O{Mt). 

The LO expression is 


Xo{xt) = ^ 


xt + 2 3xt-6 

+ 7 - TWT log Xt 


Xt-l ' {xt - 1)2 
The NLO QCD correction (ij-lsl reads, in the MS scheme 


29xt - Xt - 4xt Xt + 9xt - Xt - Xt , 
Xi{xt) = — - ,9 — • - — ^ log Xt 


3(1 - xQ2 


(1 - Xtf 


+ 


8xt + iXt + Xt — Xt 2 „ ^Xt — Xt 


2(1 


log Xt - 


(1 - Xt)2 


Li2(l - Xt) 


, . dXo g. 
+ 8xt ^— log 


dxt 




(A.2) 


(A.3) 
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where [i is the renormalisation scale. The 2-loop EW correction Xew has been 
calculated in [^, but no explicit result has been presented. Approximate formulae, 
one of which accurate to more than 0.05%, as well as a plot of the contribution 
(a/47r)AEw can however be found in [^. 

The left panel of Figure shows a plot of the scale- and scheme-independent 
quantity 


A(^) 


a(/x) sin^ 
sin^ oi{Mz) 




(A.4) 


as a function of the renormalisation scale /r, together with the la bands correspond¬ 
ing to the theoretical error in the matching of the top Yukawa coupling yt at the 
weak scale, and the experimental error on the top mass Mj. The MS couplings at 
full NNLO - 2-loop matching at the weak scale (3-loop QCD for as and yt) and 
3-loop running (4-loop QCD for ag) - as determined in have been used. The 
remaining scale dependence shown in the figure comes from higher order corrections, 
mainly from QCD, and accounts for an error of 0.004 on Xt- An additional error of 
0.002 comes from the ambiguity in the choice of the renormalisation scheme for the 
EW prefactor, as shown in [^. A comparison of the different errors contributing to 
Xt is shown in the right panel of Figure]^ The experimental error on the top quark 
pole mass Mt is by far the dominant contribution at present. 


B Expression for Pc{X) 

An approximate formula for i-(?^(A) taken from reads 

( X \ 0.5081 / tnx \ \ 1-0192 
mc(mc) \ ( as[Mz) 


P^^{X) = 0.38049 


± 0.008707 


1.30 GeV 


0.1176 


1 + ^ 




, X X 0.5276 / tux \ \ 1-8970 
mcjmc) \ ( as[Mz) \ 

1.30 Gevy V 0-1176 y 


1 + > 
(B.l) 




where 


and 


— In 


rricimc) 
1.30 GeV 


) ^a. 


Ln, = In 


asjMz) 

0.1176 


(B.2) 


Kio = 1.6624, Koi = —2.3537, kh = —1.5862, K 20 = 1.5036, kq2 = —4.3477, 

eio = -0.3537, eoi = 0.6003, en = -4.7652, 620 = 1.0253, eo 2 = 0.8866. 

(B.3) 


C More details on e'/e 


The basic one-loop functions entering rt5.6^ are given by (A.2) and 


ay 

8 


Xt-4: 

xt-1 


3xt 

{xt - 1)^ 


In Xt 


yo{xt) 


(C.l) 
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Total error: 0.6% 

1.50 

1.49 


1.48 


1.47 


100 150 200 250 300 350 

[GeV] 

Figure 8: Left: Renormalisation scale dependence of the quantity X(/x). The dashed lines 
show the uncertainty due to the error on the measured pole top mass Mt, while the dotted 
lines correspond to the theoretical error on the MS top Yukawa coupling yt due to higher 
orders in the matching at the weak scale. Right: different sources of error affecting Xt. 



^ , 1, 184 - 163x? + 259xf - 108xt 

Z„(.,)= - - In.. + i44(.. - 1)3 -+ 


+ 


324 - 38x? - 15xf + 18xi 
72{xt - 1)4 


In Xt 


(C.2) 


r:, r 2 x|(15 - 16xt+ 4x|) Xt{l8 - llxt - x"^) 

£■„(.,) = -pn., +- 

where xt = m^/M^ 


w 


The coefficients rf'\ rf''^ and rf’'^ entering (|5.7|) are given in the NDR renormal¬ 


isation scheme for p = nic and three values of as{Mz) in Table 
(1/2) 3/2 

The parameters Bq ' ' and Bg are related to the hadronic matrix elements Qq 
and Qs as follows 


(Q6(p))o — - 4 


m 


-I 2 


K 


msiy) + md{y) 


{Fk - Ft,) B^ 


( 1 / 2 ) 


(C.4) 


{Qs{y))2 = V2 


m 


K 


Fn B. 


(3/2) 


(C.5) 


rrisiy) + mdiy) _ 

It should be emphasised that the overall factor in these expressions depends on the 
normalisation of the amplitudes Ro, 2 - The matrix elements given above correspond 

On the other hand the RBC-UKQCD 


to the normalisation used in |23, 46 93 


collaboration 82,94 uses a different normalisation adopted in 91. By comparing 


(C.4) and (C.5) with Eqs. (5.10) and (5.18) of the latter paper we find that the 

have and additional factor of 4^/2. While e'/e clearly does 


91 


matrix elements in 

not depend on this difference, it is crucial to take it into account when extracting 

('3/2') 

the value oi Bg ^ from the results obtained by RBC-UKQCD collaboration. To 
this end we use Eq. (30) for A 2 in [94], adjust to our normalisation, and compare 
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as{Mz) = 0.1179 

as{Mz) = 0.1185 

as{Mz) = 0.1191 

i 

' i 

' i 

' i 

^(0) 

' i 

' i 

' i 

' i 

' i 

' i 

0 

-3.392 

15.293 

1.271 

-3.421 

15.624 

1.231 

-3.451 

15.967 

1.191 

^0 

0.655 

0.029 

0. 

0.655 

0.030 

0. 

0.655 

0.031 

0. 

>0 

0.451 

0.114 

0. 

0.449 

0.116 

0. 

0.447 

0.118 

0. 


0.406 

-0.022 

-13.435 

0.420 

-0.022 

-13.649 

0.435 

-0.023 

-13.872 

Eo 

0.229 

-1.760 

0.652 

0.228 

-1.788 

0.665 

0.226 

-1.816 

0.678 


Table 4: The coefficients rf \ rf'^ and of formula (5.1) in the NDR-MS scheme for 
three values ofas{Mz). From (83l 


to A 2 expressed in terms of (Q 8 (Ai ))2 in (C.5). This a. 
to the hadronic matrix element used in 

io>oMix 


82 


lows US to related (( 58 (m ))2 
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{Q%{h))2 = 


MS-NDR 


3\/2 


In this manner we find 
?(3/2) 




1 

6i7r 


msilj) + mdin) 


1 2 


m 


K 


a^MS-NDR 

‘( 8 , 8 )„ 


ih) 


(C.6) 


(C.7) 


The a dependence of is practically cancelled by the one of quark 

(3/2) - 

masses so that is practically ^-independent. In particular in the MS-NDR 

scheme the /r-dependence is very weak [M] . 


Using the QCD lattice value from 82|^ 


GeV) = 4.55 ± 0.27, (C.8) 

together with the light quark mass values 

ms{2 GeV) = (93.8 ± 2.4) MeV, md{2 GeV) = (4.68 ± 0.16) MeV, (C.9) 
we find 


Rf/^^(3GeV) = 0.75 ±0.05, 


B, 


(3/2) 


(me) = 0.76 ±0.05. 


(G.IO) 
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